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ABSTRACT 

We present the Hubble diagram (HD) of 66 Gamma Ray Bm^sts (GRBs) derived 
using only data from their X-ray afterglow lightcurve. To this end, we use the recently 
updated Lx - Ta correlation between the break time Ta and the X - ray luminosity Lx 
measured at Ta calibrated from a sample of Swift GRBs with lightcurves well fitted 
by the Willingale et al. (2007) model. We then investigate the use of this HD to 
constrain cosmological parameters when used alone or in combination with other data 
showing that the use of GRBs leads to constraints in agreement with previous results in 
literature. We finally argue that a larger sample of high luminosity GRBs can provide 
a valuable information in the search for the correct cosmological model. 
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1 INTRODUCTION 



That the universe is spatially flat, has a subcritical mat- 
ter content and is undergoing a phase of accelerated expan- 
sion are nowadays globally accepted ideas. Strong evidence 
for this scenario come from the anisotropy and polariza- 
tion spectra of the Cosmic Microwave Background Radia- 
tion (CMBR) (de Bernardis et al. 2000; Brown et al. 2009; 
Komatsu et al. 2010), the galaxy power spectrum (Tegmark 
et al. 2006; Percival et al. 2007) with its Baryonic Acoustic 
Oscillations (BAO) (Eisenstein et al. 2005; Percival et al. 
2010) and the Hubble diagram (HD) of Type la Supernovae 
(SNela) (Kowalski et al. 2008; Hicken et al. 2009; Kessler et 
al. 2009). What is driving this cosmic speed up and domi- 
nating the energy budget is, on the contrary, a still hotly de- 
bated question with answers running from the classical cos- 
mological constant (Carroll et al. 1992; Sahni & Starobinsky 
2000), scalar fields (Peebles & Ratra 2003; Copeland et al. 
2006) and higher order gravity theories (Capozziello & Fran- 
caviglia 2008; Sotiriou & Faraoni 2010; Nojiri & Odintsov 
2008; De FeUce & Tsujikawa 2010). In order to break the 
degeneracies among model parameters and go to the next 
level of model selection, it is mandatory to probe the back- 
ground evolution up to very high redshift so that the onset 



of the acceleration epoch and the transition to the matter 
dominated regimes may be directly investigated. 



Thanks to their enormous energy release, GRBs are vis- 
ible up to very high z, the largest one being at z = 8.2 (Sal- 
vaterra et al. 2009), hence appearing as ideal candidates for 
this task. Unfortunately, GRBs are everything but standard 
candles because their peak luminosity spans a wide range. 
There have nevertheless been many attempts to make them 
standardizeable candles resorting to the use of empirical 
correlations among distance dependent quantities and rest 
frame observables (Amati et al. 2008; Fenimore & Ramirez - 
Ruiz 2000; Norris et al. 2000; Ghirlanda et al. 2004; Liang 
& Zhang 2005). Such empirical relations allow one to infer 
the GRB rest frame luminosity or energy from an observer 
frame measured quantity so that the distance modulus can 
be estimated with an error mainly depending on the intrinsic 
scatter of the adopted correlation. Combining the estimates 
from different correlations, Schaefer (2007) first derived the 
GRBs HD for 69 objects, while Cardone et al. (2009) used 
an enlarged sample and a different calibration method to 
update the GRBs HD. Many attempts on using GRBs as 
cosmological tools have since then been performed (see, e.g., 
Firmani et al. 2006, Liang et al. 2009, Qi & Lu 2009, Izzo 
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Figure 1. Best fit caUbration parameters (a, fe, Uint) and rins of the residuals for the Lx -Ta correlation fitted to the 8 canonical GRBs 
as a function of the matter density parameter ^m for three different values of uiq, namely uiq = —1.25 (short dashed), uiq = —1.0 (solid), 
«!o = —0.75 (long dashed). We use a CPL model with h = 0.70 and set Wa = —wq to have a matter dominated era at high redshift. 



et al. 2009 and refs. therein) showing the interest in this 
application of GRBs. 

In this Letter, we rely on the Lx - Ta correlation (Dain- 
otti et al. 2008) to build the HD of 66 GRBs from their X- 
ray afterglow lightcurves observed with the Swift satellite 
(Evans et al. 2009). We then present a preliminary applica- 
tion of the derived GRBs HD showing how this dataset can 
help constrain the parameters of some simple dark energy 
models. 



2 THE Lx - Ta CORRELATION 

We use the updated version of the Lx - Ta correlation be- 
tween the luminosity Lx at the break time Ta and Ta it- 
self. We first remind the reader that Ta is defined as the 
time marking the passage from the plateau phase to the 
power -law decay in the GRBs X-ray afterglow lightcurve 
as described by the universal fitting function proposed by 
Willingale et al. (2007, hereafter W07). Using a sample of 
34 GRBs with lightcurve measured by the Swift satellite, 
Dainotti et al. (2008) first discovered that Lx and Ta are an- 
ticorrelated, as later confirmed by the semiempirical models 
of Ghisellini et al. (2009) and Yamazaki (2009). Recently, we 
have increased the GRBs sample and rederived the Lx - Ta 
correlation (Dainotti et al. 2010, herefater DIO). Introduc- 
ing the error parameter u = -v/o^^ + a'^ , DIO have also 
selected a class of high luminosity long GRBs with very well 
measured {Lx,Ta) parameters {u < 0.095) and lightcurve 
closely matching the W07 model. Referring to this class of 



objects as canonical GRBs, DIO have demonstrated that 
they define an upper envelope for the Lx - Ta correlation 
with the same slope, but a higher intercept than the one for 
the full sample. In order to avoid mixing objects possibly 
belonging to two different classes, we divide the 66 GRBs in 
two subsamples according to the value of u being smaller or 
larger than 0.095. In the first case, we select 8 GRBs which 
we will refer to as the canonical (C) sample, while the re- 
maining 58 will form the non canonical (NC) one. Adopting 
a flat ACDM model with {nM,h) = (0.278,0.699), we use 
the Bayesian method of D' Agostini (2005) to fit a linear 
relation, log Lx ~ a log [Ta/(1 + z)] + b. We thus get : 



logLx = -1.00 log 



Ta 



49.26 



.l + z 

with intrinsic scatter aint = 0.66 for the NC sample and 
T 



(1) 



logLx = -1.04 log 



l + z 



+ 50.22 



(2) 



with (Tint = 0.23 for the C sample. Note that these are the 
best fit values which have to be used in the later determina- 



* The D' Agostini method provides a well motivated approach 
to deal with the problem of fitting a linear relation when the 
uncertainties on both the {x,y) variable are comparable. Such a 
linear relation can be the outcome of a theoretical model so that 
one can expect that deviations from the underlying assumptions 
lead to the point scattering around the best fit line. This is what 
we refer to as intrinsic scatter. The D' Agostini method allows to 
take care of this term and estimate it in an unbiased way. 
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tion of the GRB distance modulus. However, the Bayesian 
approach allows to get the constraints on each one of the 
single [a, b, (Tint) by marginalizing over the other two. We 
thus find : 



-I r,n+0-2i +0-42 
-i.UU_o.21 -0.42 



. 40 Ofl+0'70 +1.41 
4y.^0_0.7o -1.37 



(Tint 



0.70 



+0.13 +0.30 



OU.iy_o.76 -1.50 



for the NC sample and 

" — J^-U4_o.22 -0.65 5 f' 

_ A .^1+0.38 +1.18 
(Tint — ".'li.o.ie -0.25 

for the C sample, where we have reported the median valueT 
and the 68 and 95% confidence ranges. In agreement with 
D09, we find that the slope is the same for the two samples, 
but the canonical GRBs are shifted to higher luminosities 
thus giving a larger zeropoint and defining an upper en- 
velope for the Lx - Ta correlation. Note that the intrinsic 
scatter for the canonical GRBs is much smaller thus leading 
to lower uncertainties on the estimated distance modulus. It 
is worth stressing that, although derived assuming a fiducial 
ACDM model, this has a negligible impact on the determi- 
nation of the distance modulus (Cardone et al. 2009). To 
strengthen this result, we have determined the best fit pa- 
rameters {a,b,aint) and the rms of the residuals for the C 
sample assuming a GPL (Chevallier & Polarski 2001; Linder 
2003) model and varying the parameters {i}M,wo,h) while 
Wa = —Wo (see later for the definition of these quantities). 
Fig. 1 shows that the best fit calibration parameters have 
a clear trend with {Qm,wo), but the end to end variation 
of (a, 6), which enter the distance modulus estimate, is less 
than 3%. This is well within the uncertainties on the (a, 6) 
coefficients so that we are confident that the choice of the 
fiducial cosmological model used in the calibration proce- 
dure has a negligible impact on the derivation of the Hubble 
diagram. However, should future data allow us to increase 
the precision on (a, 6), one should likely readdress this prob- 
lem looking for a model independent calibration. 

In order to infer the distance modulus of each GRB, 
we then simply note that Lx is related to the luminosity 
distance dhiz) as (Dainotti et al. 2008; Gardone et al. 2009) : 



Lx=47rdi(z)(l + z) 



-(2+/3) 



Fx 



(3) 



with j3 the slope of the energy spectrum (modelled as a sim- 
ple power- law) and Fx the observed flux both measured 
at the break time Ta- Having measured {Ta, l3,Fx) and in- 
ferred Lx using Eq.(l) or (2), we can then estimate the GRB 
distance modulus as : 



fi{z) = 25+5 log dL(z) = 25+- log 



47r(l + 2)-(2+/3)Fj^ 



(4) 



where dL^z) is in Mpc. The uncertainty is estimated by prop- 
agating the errors on (/3, Fx, Lx)- Note that, when comput- 
ing the error on Lx, we add in quadrature the uncertainty 
coming from Ta and Uint to take care of the intrinsic scatter. 

T Note that, because of parameters degeneracies, the marginal- 
ized likelihoods are not symmetric functions so that the median 
value may differ from the maximum likelihood one and the confi- 
dence ranges be asymmetric. See Dainotti et al. (2008) and refs. 
therein for a discussion of this issue. 
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Figure 2. GRBs Hubble diagram for the combined NC + C (black 
points) and C (red points) samples. 



When resorting to Eq.(4), we take care of the different cali- 
bration parameters for the two subsamples thus using Eq.(l) 
for the objects in the NC sample and Eq.(2) for the G sam- 
ple ones. The combined HD, shown in Fig. 2, covers the wide 
redshift range (0.033, 8.2) thus showing that, although the 
error bars at the moment are still quite large, the Lx-Ta 
correlation could allow us to probe of both the dark energy 
epoch and the matter dominated era with a single tracer. 



3 GRBS HD AS A COSMOLOGICAL TOOL 

Whatever is the tracer used, the HD is a primary tool to 
investigate the viability of a cosmological model. Indeed, 
the luminosity distance dL{z) reads : 



dLiz) 



Ho 



(1 



dz' 



(5) 



E{z') 

with E{z) = H{z)/Ho the dimensionless Hubble parameter. 
For a spatially flat universe made out of dust matter and 
dark energy with the GPL equation of state (eos) w(z) = 
Wo + Waz/{1 + z), it is: 

3WaZ^ 



E 



nM{l+zf + {^-^M){l+z)'"-'+-°+-^'> exp (-f^) (6) 



where flm is the present day matter density parameter and 
one recovers the AGDM model for (wojWa) = (—1,0). In 
order to constrain the model parameters, we use a Markov 
Ghain Monte Garlo (MGMG) algorithm to maximize the 
likelihood function C{p) oc exp [— x^(p)/2] where p is the set 
of model parameters and the expression for x^ (p) depends 
on the dataset used. Since we are here interested in testing 
the usefulness of the Lx - Ta correlation as a cosmological 
tool, we consider, as a first test, GRBs only thus setting : 



J^GRt 



x^(p) = E 



yiohs{Zi) — fJ,th(Zi,p) 



(Ti 



/i - 0.742 \^ / gjAi- 0.1356 y 



0.036 



0.0034 / 



(7) 



Here, Hots and nth are the observed and theoretically pre- 
dicted values of the distance modulus, while the sum is over 
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Table 1. Constraints on the model parameters using GRBs and priors on (h,uiM)- Columns are as follow ; 1. id of the model, 2. maximum 
likelihood parameters, 3. ,4. ,5. ,6. median and 68 and 95% CL after marginalization. Upper (lower) half of the table refers to the results 
using the NC + C (C only) sample. A sign "-" means that the parameter is fixed to its theoretical value. 



Id 



Pml 



f^A 



h 



ACDM (0.248,0.739) 

QCDM (0.248,-0.37,0.741) 

CPL (0.242,-0.34,0.02,0.751) 



n 9/IR+U-022 +0.050 
U.zy:D_Q Q24 -0.039 



0.248 



+0.026 +0.052 
-0.024 -0.042 



n 9A7+0-029 +0.066 
"■^^'-0.023 -0.043 



+0.14 +0.19 



-0.52 



_n 4S+0-11 +0-14 
^■^°-0.38 -0.96 



0.742 



+0.040 +0.068 
-0.032 -0.066 



n 7OQ+0.035 +0.067 
U. (.3M_Q Q34 _o.067 



-n fi1+l-20 +2.43 n. YAn+0-037 +0.070 
^■"^-1.39-2.21 "-'^^-0.038-0.080 



ACDM (0.246,0.743) 

QCDM (0.244,-1.26,0.745) 



0.246 



+0.027 +0.057 
-0.020 -0.037 



n 9Afi+0-'326 +0.055 
U.ztDjj (J24 -0.041 



CPL (0.245,-0.36,-2.89,0.744) 0.242+°;°^^ t°;°^^ 



-0.73 



+0.35 +0.40 
-0.72 -0.91 



n 70+0.27 +0.38 
^■'■^-0.43 -0.85 



n ni+l-44 +2.55 
"•"-'--1. 90 -2.82 



n 7/1Q+0.031 +0.061 
"■'^"'-0.039 -0.067 

n 7/iq+0.036 +0.074 
"•'^'^-0.036 -0.069 



n 7/1S+0-034 +0.077 
"•'^°-0.036 -0.080 



the Ngrb GRBs in the sample. The last two terms are Gaus- 
sian priors on h and lom = 0,Mh^ and are included in order 
to help break the degeneracies among model parameters. To 
this aim, we have resorted to the results of the SHOES col- 
laboration (Riess et al. 2009) and the WMAP7 constraints 
(Komatsu et al. 2010), respectively, to set the numbers used 
in Eq.(7). We consider three different cosmological models, 
namely the ACDM one, quiessence (QCDM) obtained by 
setting Wa = but wo free to vary, and the full CPL model 
leaving both (wo, Wa) unspecified. The results of the MCMC 
analysis are summarized in Table 1 for the three models con- 
sidered and two GRBs samples, namely the full NC -\- C set 
or the C subsample only. Note that this test allows us to 
investigate whether a sample made of only high precision 
GRBs is better suited to constrain cosmological parame- 
ters. As a general remark, we note that the constraints are 
in good agreement with previous ones in literature (Kowalski 
et al. 2008; Hicken et al. 2009; Komatsu et al. 2010; Percival 
et al. 2010) thus showing that using GRBs does not intro- 
duce any bias in the estimate of cosmological parameters. 
While the constraints on {Q.m, h) are quite narrow, those on 

the eos parameters are quite weak+ when deviations from 
the ACDM model are allowed. This is not surprising since 
the GRBs HDs mainly probe the matter dominated high - z 
regime. As such, the details of the eos have only a minor 
impact on the background evolution and the HD. We nev- 
ertheless note that, although with a large error, the fit to 
the CPL model points towards Wa = 0, i.e. a constant eos. 
Considering that wo = — 1 is well within the 68% CL range, 
one can argue that the ACDM model is still preferred by a 
GRBs only fit. Such a conclusion is also strengthened not- 
ing that the reduced x^ (i-e., X^ /{^grb — rip) with rip the 
number of parameters) for the maximum likelihood model 
is the same for the three cases considered so that increasing 
rip is not statistically supported. 

■f Note that the 68% and 95% CL on uiq in Table 1 are essentially 
driven by the assumption wq < ~l/3 we have imposed in order 
to have p + 3p < for the dark energy fluid. 



It is worth comparing the results from the NC + C and 
C samples. The canonical GRBs cover a narrower redshift 
range (0.125 < z < 2.612 (with a median Zmed ~ 1.26) than 
the NC sample (with 0.033 < z < 8.2 and z,ried = 2.31), but 
with significantly less uncertainties on the distance modulus 
estimates. In a sense, one can say that comparing the results 
from the two samples helps understanding which strategy 
(increasing statistics or ameliorating precision) should be 
chosen to improve the constraint on cosmological parame- 
ters by GRBs. Comparing the values in the upper and lower 
half of Table 1, one sees that the percentage error on each 
parameter is comparable using either the NC -I- C or the C 
sample, while the constraints on the parameters (in particu- 
lar, flM and h) are in remarkable good agreement. One can 
therefore argue that the reduction in the number of GRBs 
from the NC -I- C to the C sample has been compensated by 
the increased accuracy, but a more extensive analysis have 
to be carried out before a definitive answer could be drawn. 

As a next step, we combine the GRBs HDs with other 
data redefining £{p) as : 



^(P) 



exp(-XGfli3/2) 



(2^)A^GflB/2|CG«s|l/2 




1 


1 / ft - hobs \ ''' 




y'2nal 




exp(-xlAo/2) 




(27r)-^BAo/2|Cs^o|i/2 


1 


1 ^TZ-TZobsY' 

2 I aij ) 


V27ra| 


exp (-xlf /2) 





(27r)A^«/2jC^|l/2 



(8) 



Here, the first two terms are the same as in Eq.(7) 
with Cgrb the GRBs diagonal covariance matrix and 
(hobsji^h) ~ (0.742,0.036). The third term takes into ac- 
count the constraints on d^ — rs{zd)/Dv{z) with rs{zd) the 
comoving sound horizon at the drag redshift Zd (which we 
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Table 2. Same as Table 1 but using GRBs + SHOES + BAO + 11 + H{z). 



Id 



Pml 



wo 



ACDM 

QCDM 

CPL 



(0.265,0.709) 
(0.271,-0.98,0.707) 



n 9R'5+0.026 +0.054 
U.ZDO_Q^Q2i _().04i 



0.267 



+0.035 +0.072 
-0.032 -0.061 



(0.262, -1.32, 1.01, 0.726) 0.273rg ' 



+0.036 +0.069 
036 -0.064 



-1.00 



+0.14 +0.29 
-0.15 -0.37 



-1.01 



+0.47 +0.64 
-0.33 -0.54 



n 700+0-015 +0.034 
U.(UJ_Q Qj^4 _o.030 



n 7in+0.025 +0.057 
"■'-'-^-0.025 -0.051 



-n n9+0'*** +1.57 r, 7r.p+0.033 +0.062 
^■'^^-1.45 -2.44 ^-'^"-0.032 -0.053 



n ofiQ+0.022 +0.049 
U.^DCi_g Q24 -0.048 



ACDM (0.262,0.710) 

QCDM (0.244,-1.26,0.745) 

CPL (0.249,-1.50,1.40,0.742) 0.256to:o32 



n 9Afi+0-026 +0.055 
"■^^"-0.024 -0.041 



+0.037 +0.085 
0.060 



-0.73 



+0.35 +0.40 



-1.06 



+0.43 +0.64 
-0.38 -0.54 



p, „,n+0-01f> +0-028 

'^■' 1*^-0.016 -0.031 

n 740+0.037 +0.074 
"•'^•^-0.035 -0.069 

-n 09+1-06 +1-54 r, 7-17+0.033 +0.059 
'-'■'^^-1.35 -2.49 ^■'-'-'-0.034 -0.064 



fix to be rs{zd) = 152.6 Mpc from WMAP7) and the volume 
distance is defined as (Eisenstein et al. 2005) : 



Dv{z) = 



H{z) 



Dl{z) 
1 + z 



1/3 



(9) 



The values of d^ ai z — 0.20 and z — 0.35 have been 
estimated by Percival et al. (2010) using the SDSS DR7 



galaxy sample so that we define Xbao ~ D 



Cg\oC with 

D^ = (dgl'l - d'o\,dtk - d'o^s) and Cbao the BAO covari- 
ance matrix. The next term refers to the shift parameter 
(Bond et al. 1997; Efstathiou & Bond 1999) : 



7^: 



dz' 



(10) 



with Zi, — 1090.10 the redshift of the last scattering 
surface. We follow again WMAP7 setting {TZobs,(yiz) ~ 
(1.725,0.019). While all these quantities (except for the 
Gaussian prior on h) mainly involve the integrated E{z), the 
last term refers to the actual measurements oi H{z) from the 
differential age of passively evolving elliptical galaxies. We 
then use the data collected by Stern et al. (2010) giving the 
values of the Hubble parameter for A/if = 11 different points 
over the redshift range 0.10 < z < 1.75 with a diagonal co- 
variance matrix. The results obtained fitting this combined 
dataset are summarized in Table 2. Not surprisingly, we find 
that the ACDM model is still statistically favoured by the 
combined dataset. The constraints on {i}M,h) for this case 
are almost unchanged by the addition of the other data. 
This is expected since the data we have added help break 
the degeneracy between the matter content and the eos pa- 
rameters. Since in the ACDM case the eos is set from the 
beginning, the priors on {h,ujM) we have used to get Table 1 
are essentially equivalent to the full dataset in Table 2. The 
situation is different for the QCDM and CPL models since 
the eos is allowed to vary. The combined dataset is now able 
to constrain (wo,Wa) thanks to the H{z) and BAO data 
which probe the redshift range where dark energy drives 
the background evolution thus complementing GRBs. In or- 
der to better show the impact of GRBs, we have repeated 
the fit excluding their HD. Considering, as a test case, the 



CPL model, we find for the maximum likelihood param- 
eters {Q,M,wo,Wa,h) = (0.258,-1.52,1.51,0.739) close to 
that obtained using the GRBs. The marginalized constraints 
now read : 



O _ n oRQ+O-033 +0.062 
iLM — U.z;Dd_Q 034 -0.061 

,,, _ _n 99+1-30 +1-82 



Wo 



-1.00 



+0.35 +0.56 
0.42 -0.60 



h — n 714+0.031 +0.056 
n — U.I 14_o.029 -0.053 



Comparing these values to those in Table 2 shows that 
adding GRBs does not significantly narrow the parameters 
confidence ranges. This is actually expected since the present 
GRBs datasets are affected by large errors (the NC sample) 
or few statistics (the C sample). However, it is worth noting 
that GRBs help push the constraints on Wa towards thus 
suggesting that the inclusion of a large sample of canonical 
{u < 0.095) GRBs may strengthen the case for a constant 
eos dark energy model. 

Up to now, we have not used SNela since we have been 
mainly interested in investigating how GRBs (and the Lx - 
Ta correlation) could be used as cosmological tools. Having 
demonstrated that the inclusion of GRBs does not bias the 
search for cosmological parameters, we finally add SNela to 
the above combined dataset using the Constitution sample 
(Hickcn et al. 2009) comprising 397 objects over the range 
0.015 < z < 1.551. For the CPL model and the NC-^C 
GRBs, we get (Om, wo, Wa, ft) = (0.279,-1.00,0.16,0.708) 
as maximum likelihood parameters, while the marginalized 
constraints on the single quantities are as follows : 



O — n 078+0-021 +0.054 
liM — U./<8_o.022 -0.041 



1 nn+0.15 +0.32 
Wo = -1.00_o.l7 -0.30 . 



0.09 



+0.68 +1.16 
-0.62 -1.69 ) 



h = 0.703 



+0.019 +0.035 
-0.015 -0.032 



Not surprisingly, the maximum likelihood model is close to 
the ACDM {wo, Wa) ~ (—1, 0) case, while the constraints on 
the parameters are in full agreement with those reported in 
Tables 1 and 2 being narrower thanks to the SNela data. In 
particular, the eos parameters are now better constrained 
since the SNela HD mainly probes the low redshift dark 
energy dominated universe. We have not tried to repeat this 
analysis using the C sample since, due to their high number 
and good accuracy, SNela actually dominate the fit as can 
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be checked by fitting SNela only witliout tlic inclusion of 
GRBs or tlie otlier data. 



4 CONCLUSIONS 

A ground zero test of every cosmological model is its abil- 
ity to fit the observed Hubble diagram. Such a test tells 
us whether the model is able to give rise to the observed 
background evolution and select the region of the parameter 
space leading to the correct sequence of accelerating and de- 
celerating expansion. Although still being the primary tracer 
of the Hubble diagram, SNela may trace cosmic evolution 
only up to 2 ~ 2 so that it is mandatory to look for a differ- 
ent class of astronomical objects to go beyond this limit and 
investigate the matter dominated era. GRBs are considered 
the ideal candidates for this role thus motivating the hunt for 
a calibration method to make them standardizeable candles. 
We have presented here the use of the Lx - Ta correlation as 
a valid method to infer the GRBs HD up to z = 8.2 using 
data coming from the X-ray lightcurve. It is worth stress- 
ing that this is the only empirical law relating quantities 
measured from the afterglow lightcurve rather than being 
related to the prompt emission quantities. Moreover, differ- 
ently from what has yet been done in the past (Schaefer 
2007; Cardone et al. 2009), the HD for the NC -|- C sample 
is the only GRBs HD based on the use of a single correlation 
and containing a statistically meaningful number of objects. 
The use of the Lx - Ta correlation then avoids the need of 
combining different correlations to increase the number of 
GRBs with a known distance modulus. Each correlation is 
affected by its own possible systematics and characterized 
by different intrinsic scatter so that combining all of them 
in a single HD can introduce unexpected features and hence 
bias the constraints on the cosmological parameters. 

The intrinsic scatter of the Lx - Ta correlation may be 
significantly reduced if one considers only its upper envelope 
defined by the canonical GRBs. On the other hand, using 
the full 66 GRBs, we have been able to get constraints on the 
matter content and the present day Hubble constant compa- 
rable to those yet available in literature. Such a preliminary 
investigation has convincingly shown that the Lx - Ta cor- 
relation may be used to construct a GRBs HD which does 
not introduce any bias in the search for cosmological pa- 
rameters. It is worth stressing that the same results have 
been obtained using the canonical GRBs despite the fact 
that they represent just ~ 12% of the full sample. We there- 
fore argue that an observational effort dedicated to look for 
canonical GRBs may turn them from ideal candidates to ac- 
tual tools for efficiently investigating the dark energy puzzle. 
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